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Discontinuities in elastic constants are detected at the superconducting transition of layered or- 
ganic conductors k-(BEDT-TTF)2X by longitudinal and transverse ultrasonic velocity measure- 
ments. Symmetry arguments show that discontinuities in shear elastic constants can be explained 
in the orthorhombic compound only if the superconducting order parameter has a mixed character 
that can be of two types, either Aig + Big or B2g + Bsg in the classification of irreducible represen- 
tations of the orthorhombic point group D2h- Consistency with other measurements suggests that 
the Aig + Big (dxy + d^^x+y)) possibility is realized. Such clear symmetry-imposed signatures of 
mixed order parameters have not been observed in other superconducting compounds. 

PACS numbers: 74.70.Kn,74.25.Ld,74.20.Pp 



Unconventional, non s-wave, superconductors in solids 
seem ubiquitously associated with strong electronic cor- 
relations. This is the case in a wide variety of compounds 
that include heavy fermions, ruthenates, cuprates as well 
as quasi-two-dimensional half-filled organic charge trans- 
fer salts K-(ET)2X (ET = BEDT-TTF) [l|. In most cases 
gaps with nodes are observed, but the exact symmetry 
of the unconventional superconducting order parameter 
is uncontroversial only in the cuprates. 

In this letter, we focus on the layered organics that 
exhibit antiferromagnetism and Mott insulating behav- 
ior, as the cuprates, and establish the two-component 
nature of the singlet order parameter in the orthorhom- 
bic compound K-(ET)2Cu[N(CN)2]Br. Previous stud- 
ies suggest d-wave pairing with nodes, although s-wave 
symmetry is sometimes seen. Measurements sensitive 
to the fc-space dispersion, such as scanning tunnel- 
ing spectroscopy and thermal conductivity favor 
dxy symmetry, namely nodes along the nearest neigh- 
bor bonds (or equivalenty, between the orthorhombic 
axes). Moreover, theoretical calculations based either on 
spin-fluctuation mediated superconductivity 0, H, H, 0] 
or on quantum cluster methods Q and variational 
approaches [l3| for the Hubbard model, support the 
anisotropic d-wave picture with a prevailing dj.2„j,2 sym- 
metry. Nevertheless, none of these calculations has con- 
sidered interlayer hopping, which, as we will show, is nec- 
essary to explain the experimental data that we present. 

The ultrasonic probe is extremely sensitive to gap 
anisotropics as the attenuation and velocity depend on 
both the direction of wave propagation and the direc- 
tion of polarization. Attenuation experiments on UPta 
[ill [l^ and of Sr2Ru04 [13,] perfectly illustrate how the 
unconventional gap structure can be unraveled by such a 
versatile technique. In organic charge transfer salts how- 



ever, attenuation experiments are hampered by the small 
size and the shape of single crystals. Nevertheless, one 
experiment was successful for the K-(ET)2Cu[N(CN)2]Br 
compound jl^] , but the interpretation of the results was 
complicated by a phase separation occurring even in 
highly ordered samples. Notwithstanding these difficul- 
ties, ultrasound velocity can also be used to obtain in- 
sights into the nature of the superconducting (SC) state 
in layered organics. Lattice anomalies [lB| and elastic 
constant changes [IB, [l3| have been identified, but no 
consistent effort has been yet dedicated to identify the 
SC order symmetry. 

We report anomalies observed at the SC transi- 
tion temperature Tc on three elastic constants of 
monoclinic k-(ET)2Cu(NCS)2 and of orthorhombic k- 
(ET)2Cu[N(CN)]2Br. Even though these compounds be- 
long to different point groups, we expect similarities in 
the SC order parameters because of their nearly iden- 
tical electronic properties. To understand discontinu- 
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Since we focus on symmetry prop- 
erties, a Ginzburg-Landau (GL) approach will suffice 

We use an acoustic interferometer [ij] to measure rel- 
ative changes in velocity AV/V that allow us to extract 
the corresponding relative variations in the elastic con- 
stants C through AC/C = 2AV/V. The k-(ET)2X crys- 
tals grow as platelets containing the highly conducting 
planes whose normal is oriented along a* for monoclinic 
k-(ET)2Cu(NCS)2 and along b for the orthorhombic k- 
(ET)2Cu[N(CN)2]Br. Thus, ultrasonic plane waves can 
be propagated only along these normal directions. Pure 
longitudinal and transverse waves cannot be propagated 
along the a* axis of the monoclinic structure so, strictly 
speaking, it is not possible to measure the Cij 's indi- 
vidually as it is the case for the orthorhombic material 
[2^. However, given the layered structure and the a axis 
orientation of about 110° instead of 90° from the plane, 
we neglect, as a first approximation, the off-diagonal ele- 
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C55 (5) 


C66 (a) 


Transverse 


Cee (b) 


C44 (5) 



TABLE L Elastic constants dj with the appropriate polari- 
sation of the ultrasonic waves for two k-(ET)2X compounds. 

ments of the Cij matrix that differentiate the monoclinic 
structure from the orthorhombic one. This simpHfies the 
data treatment without affecting the conclusions. With 
this approximation the measured C^j's are given in Ta- 
ble |T] for each crystal. 
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FIG. 1: (Color online) Longitudinal waves propagating along 
the a* axis in k-(ET)2Cu(NCS)2: (A) AV/V data at 166 MHz 
for H = and 12 Tesla ; (B) ACn/Cn at three frequencies. 

The «;-(ET)2Cu(NCS)2 crystal will be considered as 
our reference compound since it is located far enough 
from the Mott transition line on the high pressure side of 
the P-T diagram with no indication of a phase separation. 
To extract the elastic change caused by the onset of su- 
perconductivity, we applied a magnetic field perpendicu- 
lar to the highly conducting plane to quench the SC state. 
We show in Fig. [T]A the temperature dependence of the 
relative change of the longitudinal velocity below 20 K at 
166 MHz. In zero magnetic field a negative discontinuity 
of the velocity is obtained at Tc = 9.5 K; the anomaly is 
completely quenched in a field of 12 Tesla leaving only 
a monotonous decrease of the velocity as the tempera- 
ture increases. We notice the absence of magnetic field 
effects above 12 K, an observation that excludes, con- 



trary to the K-(ET)2Cu[N(CN)2]Br compound [11], the 
presence of a coexisting phase in this temperature range. 
The difference between these two curves yields the rela- 
tive variation of the compressional constant Cn shown 
in Fig. [Tj3 at different frequencies. As expected, no fre- 
quency dependence is observed: the onset of the SC phase 
yields a negative discontinuity at Tc that extends over 
a few degrees due to important SC fluctuations above 
and below the superconducting temperature defined as 
the maximum slope. At lower temperatures ACn/Cn 
is practically constant. A similar procedure was used 
for the two transverse acoustic modes yielding, over the 
same temperature range, AC55/C55 and ACee/C'ee- The 
three relative elastic constant variations are compared in 
Fig. [21 While a negative discontinuity is expected on Cn 
[18| . the appearance of a discontinuity on the shear con- 
stant C55 is unusual. The amplitude of the discontinuity 
is larger than that of Cn by approximately a factor two, 
excluding the simple explanation of mode mixing for a 
quasi-transverse wave. These discontinuities are larger 
than in other non conventional superconductors [13, [3| 
by two to three orders of magnitude. No discontinuity is 
observed for ACge/C'ee; only a small change of slope is 
obtained at Tc- 



= 9.5 K 
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FIG. 2: (Color online) Temperature dependence of Adj/dj 
for the K-(ET)2Cu(NCS)2 compound. The vertical dashed 
line indicates the SC critical temperature. 

In the highly ordered K-(ET)2Cu[N(CN)2]Br com- 
pound, the coexistence of antiferromagnetic (AF) and 
SC phases complicates the analysis of the SC state [l4| . 
Moreover, a higher magnetic field is needed to quench 
the SC state. We present in Fig. [3] the Adj/Cij ob- 
tained by substracting the zero and 16 Tesla curves. We 
notice that magnetic field effects are observed in the nor- 
mal state up to 20 K on C22 and Cge (ACy /Cy is not 
zero). The temperature dependence below Tc = 11.9 
K is also not monotonous and the SC fiuctuations ap- 
pear on a wider temperature range above Tc- Notwith- 
standing these differences, the comparison with the k- 
(ET)2Cu(NCS)2 data (see Fig. [2]) at Tc is remarkable: 
we still observe a negative discontinuity on C22 (Cn), a 
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larger one on Cge (C'55) and only a change of slope on C44 
(Cge)- These observations clearly establish the similarity 
of the couplings between the SC order parameter and the 
elastic strains, although the crystal symmetry groups dif- 
fer because of the tilting of the axis normal to the planes. 
Moreover, they confirm that the negative dicontinuity on 
AC55/C55 for the monoclinic compound is intrinsic and 
that it cannot be attributed to mode mixing. 
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FIG. 3: (Color online) Temperature dependence of Adj/dj 
for the K-(ET)2Cu[N(CN)2]Br compound obtained with a 16 
Tesla magnetic field. The dashed line indicates the SC critical 
temperature. 

Experiment has established that the layered organics 
are singlet superconductors [l| . In the simplest GL model 
then, discontinuities in elastic constants at the supercon- 
ducting transition are easily explained through the free 
energy functional 

F = a\v\^ + geM^ + ^\tj\^ + ^ e.e, (1) 

where 77 is the order parameter, 6 is a constant, £i is the 
strain, Cij the matrix of elastic constants, while a is pro- 
portional to (T — Tc). If one of the strains is coupled lin- 
early through the constant g to the order parameter, the 
minimization with respect to rj shows that at the tran- 
sition a negative discontinuity appears on the effective 
elastic constant = d'^F/de-f. Such a linear couphng 
to is possible only if the strain is invariant under all 
the operations of the point group because jTyp is. Higher 
order coupling terms in the free energy would only lead 
to the change of slope or curvature observed below Tc for 
all Cii, and these are not considered here. 

Table |lT] shows a simplified character table for the ir- 
reducible representations of the monoclinic C2/1 group of 
K-(ET)2Cu(NCS)2, along with the transformation prop- 
erties of the strains and examples of basis functions for 
the order parameter. Note that the x and y axis are 
not perpendicular. They lie along the atomic bonds, 
which are along the diagonal formed by the b and c axes. 
Since, according to Table[TTl ei and £5 are invariant under 



irrep 


E C\ 


Basis functions 


Strains 


A, 


1 1 


s, xy, {x -f y)z 


ei, £2, £3, £5 


B, 


1 -1 




£4, £6 



TABLE II: Simplified character table, basis functions and 
transformation properties of the strains. The monoclinic Czh 
group for K-(ET)2Cu(NCS)2 has the character table of C2 ®i, 
but inversion i always has character -1-1 for singlets so the ta- 
ble of C2 shown above suffices. The names of the irreducible 
representations are those of C2h- The last column shows the 
transformation properties of the strains and the next to last 
column examples of basis functions for the order parameter. 
The a axis is tilted towards c (equivalently x + y) axis in the 
layers. 



irrep 


rp y^a y^b /^c 

-C/ O2 O2 O2 


Basis fcts 


Strains 


Big 
B2g 

Bag 


1111 
1-1-1 1 

I- 1 1-1 

II- 1-1 


s, xy 

{x + y)z 

2 2 
X ~y 

(x - y)z 


£1, £2, £3 

£6 

£5 

£4 



TABLE III: Simplified character table, basis functions and 
transformation properties of the strains for the orthorhombic 
-D2/1 = -D2 ® i group appropriate for fi:-(ET)2Cu[N(CN)2]Br. 
The b axis is perpendicular to the layers and the x + y axis 
is along a. 

the symmetry operations of the group, the correspond- 
ing elastic constants can couple linearly to \ri\'^, leading 
to negative discontinuities. However, Eq is not invariant 
so there is no discontinuity at Tc- This explains the ob- 
servations for K-(ET)2Cu(NCS)2 and it does not impose 
any constraint on the symmetry of the order parameter. 

In the orthorhombic K-(ET)2Cu[N(CN)2]Br, because 
of the different conventions, the role of €5 in the mono- 
clinic case is played by eg . The simplified character table 
mil for the D2h group shows that the shear strain eg is 
not invariant under the operations of the group. Hence, 
the Cee negative discontinuity at Tc cannot be explained 
with the simplest model Eq.([T]). One must introduce an 
order parameter with two orthonormal basis functions 
with respective complex coefficients 771 and 772. Let us 
first neglect the strain terms and consider the most gen- 
eral free energy functional that is invariant under the 
|)CTnt group and phase changes of the order parameter 

F, = a,\m\' + ^-^\vi\' + a2\m\' + j\m\' 

+ \m\^\V2\\j + Scos{2A0)). (2) 

In this expression, 7 and 6 are constants and A6 is the 
phase difference between the two components of the or- 
der parameter. If S is positive, this free energy will be 
minimized by A9 = ±7r/2, while if S is negative A9 = 
or TT will be the minimum. The case A9 = ±7r/2 cor- 
responds to a complex order parameter, hence it breaks 
time reversal symmetry. 
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To explain the discontinuity in the transverse elastic 
constant, the coupling free energy 



Fr,e = .ge6hi||?72|cos(A6') 



(3) 



must be allowed by symmetry. Also, cos(A6') should not 
vanish, thus removing the possibility of a time-reversal 
symmetry-breaking state. Since eq transforms according 
to the Big representation, there are only two possibilities. 
Either one of the 77 is invariant {Aig) and the other one 
transforms as Big or one of the components transforms 
like B2g and the other one like Bj,g. This can be checked 
by showing that the product of the characters in Table 
mil is unity for all group operations applied to F^^. Note 
that both of the above possibilities for rji and 772 forbid 
a linear coupling to £4 since the latter transforms like 
B^g. This explains the absence of a discontinuity in the 
corresponding elastic constant. 

Since both scanning tunneling spectroscopy [2] and 
thermal conductivity [3] suggest nodes along the x and 
y axis, this forces us to choose an order parameter that 
has a mixed Aig + Big character, namely d^y + d^^x+y)- 
The nodeless s case has the same symmetry as d^y so 
more generally it should be included but it suffices that 
its amplitude be smaller than that of dxy for the nodes 
of s -|- dxy to survive. They are just shifted from their 
position in the dxy case. The d^i^x+y) component does 
not remove the nodes in the plane, but it clearly breaks 
mirror symmetry about the planes. 

On general grounds, the free energy Eq.Q predicts 
two different T^'s since there is no a priori reason why 
ai and 02 should vanish at the same T. That is dif- 
ferent from the case of Sr2Ru04 where the two compo- 
nents of the order parameter necessary to explain the 
data belong to a single two-dimensional representation 
i?2u of the point group [12]. Although the present 
lattice is nearly triangular, the two components of the 
order parameter that we found do not coalesce into a 
sing le two-dimensional representation of the ZJg/i group 
[30]. Nevertheless, the mixed Aig + Big representation 
for the orthorhombic crystal does coalesce into the one- 



dimensional Ag representation of its monoclinic cousin, 
leading to a single Tc in that case. Hence, we do not 
expect a large difference between the two transition tem- 
peratures of the orthorhombic crystal. Our experimental 
data in Fig. [3] show a rather broad transition with an 
extended region of SC fluctuations that could mask this 
difference between the two transitions. 

The presence of a dz(x+y) component to the order pa- 
rameter suggests that interlayer hopping is an important 
variable in the problem. The value of this parameter 
has been estimated from angle-dependent magnetoresis- 
tance oscillations [3l[. Thermal expansion data 32] also 
disclose a striking anisotropy and dependence of on 
interlayer effects [33j| that are unlikely to be captured by 
a 2D purely electronic model. 

In summary, symmetry and the observed discontinu- 
ities at Tc in the ultrasonic velocity data for two com- 
pounds of the layered k-(ET)2X organic superconduc- 
tors demonstrate that the order parameter must have 
at least two components in the orthorhombic compound 
K-(ET)2CujN(CN)2]Br. Consistency with other exper- 
iments selects Aig -\- Big (equivalently dxy -f dz(x+y))- 
The two components coalesce into a one-dimensional ir- 
reducible representation Ag in the monoclinic compound 
K-(ET)2Cu(NCS)2. Nodes are not symmetry imposed 
but are symmetry allowed and are likely to occur in elec- 
tronic pairing mechanisms. The d^i^x+y) component of 
the order parameter suggests that further studies of in- 
terlayer coupling are called for. 
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